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SPECTRAL AND INTERMITTENCY PROPERTIES OF RELATIVISTIC TURBULENCE 
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ABSTRACT 

High-resolution numerical simulations are utilized to examine isotropic turbulence in a compressible 
fluid when long wavelength velocity fluctuations approach light speed. Spectral analysis reveals an 
inertial sub-range of relativistic motions with a broadly 5/3 index. The use of generalized Lorentz- 
covariant structure functions based on the four-velocity is proposed. These structure functions extend 
the She-Leveque model for intermittency into the relativistic regime. 

Subject headings: hydrodynamics — methods: numerical — gamma-ray burst: general — turbulence 



1. INTRODUCTION 

A rapidly accumulating body of astronomical obser- 
vations indicates the presence of relativistic turbulence 
in a variety of astrophysical systems. Cosmological 
gamma-ray bursts (GRBs) accelerate relativistic jets 
with energies exceeding the kinetic energy of a supernova 
(10 51-52 ergs) to Lorentz factors 7 > 100 and contain in- 
ternal fluctuations with ^7 ~ 2 (Piran 2004), while jets 
from active galactic nuclei and micro-quasars are found 
to have 7 > 10 and 7 < 2 respectively. Recent dis- 
coveries indicate that a broad class of stellar explosions, 
including those asso ciated with X-ray transients and low- 
luminosity GRBs (jKulkarni et all Il998t iStarling et"aT] 
l20H . and even some supernovae lacking high-e nergy 
emission (jSoderberg et alJl2010t iParagi et aT]l2010t l , also 
accelerate significant amounts of relativistic material. 
These relativistic flows are highly susceptible to turbu- 
lence due to the extremely high inferred Reynolds num- 
bers and the abundant presence of shear due to observed 
flow asymmetries. 

In recent models for fluctuating emission from GRBs, 
the relativistic turbulence is proposed to be directly ob- 
servable, as relativistic eddies beam their radiation into 
narrow cones which sweep by the observers' line of sight 
(N aravan fe Kumad [20091 lLazar et al.1 12001 . In addi- 
tion, the central engine for short-duration GRBs likely 
involves the turbulent merging of neutro n stars which are 
know n sources of gravitational radiation (|Hulse fc Tavlorl 
119751) . Knowledg e of the properties of turb ulence present 
in these sources (jMelatos fe Peralta|[20i0l ) may aid with 
detection and interpretation of gravitational waves from 
these sources when interferometers such as advanced 
LIGO become operable in the near future. Relativistic 
turbulence is also likely generated by phase transitions 
in the early universe and may be produced in heavy ion 
colliders and laser experiments. 

The fundamental properties of kinematically rela- 
tivistic hydrodynamic turbulence have not, until now, 
been investigated. Turbulence in full y relativistic 
MHD has been studied for decaying dZhang et all 
[20091 llnoue et all I20TH IBeckwith fe Stone) 1201 II) and 
driven (|Zrake fe MacFadvenl 120111 12012ft flows. Also, 
simulations ignoring th e uid res t mass have been 
done in the magnetized jChol 120051 ) and unmagnetized 
(jRadice fe Rezzollal[2012l ) cases. Fundamental questions 



relating to the non-linear behavior of the SRHD system 
remain to be answered. For example, it has not pre- 
viously been known whether relativistic velocity excita- 
tions survive to small scales or are suppressed near the 
outer scales of a turbulent flow. If relativistic motions 
persist deeply into the cascade, then universal behavior 
may be expected. However, relativistic statistical mea- 
sures are required to characterize its scaling and inter- 
mittency properties. For example, the power spectrum 
of three- velocity cannot obey a power law above the scale 
€ 7 where the velocity differences asymptotically approach 
light speed. Traditional two-point structure functions 
also break down at scales I > £ 7 . Statistical diagnostics 
built on the fluid four- velocity, u = (7c, 7V) , are thus 
needed. 

In this Letter, we present high-resolution simulations 
of isotropic relativistic turbulence and introduce Lorentz 
covariant scaling diagnostics to analyze its basic proper- 
ties. Questions to be addressed include (1) Does rela- 
tivistic turbulence display universal behavior? (2) Can 
its intermittency be described in terms of known phe- 
nomenological models? (3) How deeply into the cascade 
do velocity fluctuations remain relativistic? 

In order to address these questions, numerical integra- 
tions of the special relativistic hydrodynamics (SRHD) 
equations were conducted on the Pleiades supercomputer 
at the Ames research center on lattice sizes up to 2048 3 , 
utilizing more than 2,000,000 CPU-hours. The SRHD 
equations are cast in conservation form for the parti- 
cle number \7 t ,N 1 ' = and energy-momentum tensor 
V^T^ = and closed by the adiabatic equation of 
state p = pe(T — 1) where T — 4/3, e is the specific 
internal energy, p is the co-moving mass density and 
T» v = phu^u" +pg liU , and h = c 2 + e+p/p. The source 
term = pa^ — p(e/eo) includes injection of energy 
and momentum at the large scales and the subtraction 
of internal energy (with parameter eo = 0.25) to per- 
mit stationary evolution. Vortical modes at k/2-n < 3 
are forced by the four-acceleration field a M = which 
smoothly deco rrelates over a large- e ddy t urnover time, 
as described in iZrake fc MacFadvenl ([2012! ) . The numer- 
ical scheme employed is the 5th-order accurate weighted 
essentially non-oscillatory (WENO) conservative finite- 
differen ce scheme which has been e xtensively tested for 
SRHD (|Zhang fc MacFadvenl [20060 and is implemented 



2 



J. Zrake and A. I. MacFadyen 




Figure 1. Slices of the co-moving density (left) and four- velocity component out of the page (right), indicating the structure of isotropic 
relativistic turbulence at a resolution of 2048 3 . Density is scaled to the 1/2 power to elucidate structures in the evacuated regions. The 
grayscale is from minimum (darkest) to maximum (brightest). 



in the Mara code (jZrake k. MacFadvenl 120121) wi th im- 
proved smoothness indicators ( Shen fc Zhal l2010h . The 
simulations take place on the periodic cube of size L, 
and were evolved at stationarity for at least 5 large eddy 
turnover times, with the exception of the 2048 3 run which 
was evolved for 1/2 of a turnover time. Each model at 
resolution N 3 was initialized from stationary turbulence 
at resolution (N/2) 3 , after which a full turnover time is 
required to re-establish stationarity. Table Q] summarizes 
the simulation parameters for the 2048 3 run. 

2. RESULTS 

Images of the fully developed turbulent fields are shown 
in Fig. [TJ The co-moving density (left panel) con- 
tains fluctuations above and below the mean by fac- 
tors of 16. The highest density structures are sheet-like, 
but not fully two-dimensional and trace the front be- 
tween clouds of fluid colliding with relativistic velocities. 
Shocks and contact discontinuities are evident where the 
density field changes abruptly. In compressible turbu- 
lence, these shocks indicate the g eometry of the mos t 
aggressively dissipating structures (|Kritsuk et al.ll2Q07l) . 
These shocks are more intense in a relativistic gas be- 
cause their density contrasts may be arbitrarily large, 
a = j^ji + jAj-e/c 2 whereas in the non- relativistic limit 
the density ratio is bounded by ^zj- The transverse four- 
velocity component (right panel) reveals a highly struc- 
tured hierarchy of shearing motions down to the smallest 
scales. It also exhibits remarkable coherence over some 
regions, where large clouds of fluid are in relativistic mo- 
tion into (or out of) the page. Fluid elements separated 
by 1/100, 1/10 and 1/5 of the image are, on average, 
in relative motion with four- velocities of 1/2, 1 and 2 
respectively. 

Turbulent flows at sufficiently large Reynolds number 
are known to contain a range of scales over which the 



Table 1 

Simulation parameters 



Parameter 


Symbol 


Value 


Box size 


L 




Driving scale scale 




L/2 


Relativistic scale 




0.16L 


Mean pairwise-rclative Lorentz factor 


TreZ 


3.10 


Mean lab-frame Lorentz factor 


7 


1.67 


Max lab-frame Lorentz factor 


')max 


5.73 


Mean relativistic Mach number 


M 


2.67 


Max relativistic Mach number 


J^max 


22.9 


Mean pressure to mean density ratio 


p/pc 2 


0.41 


Lower 1% density quantile 




0.027p 


Upper 1% density quantile 




4.89p 



Note. — Simulation parameters are taken from the 2048 
run. The relativistic Mach number is defined to be JVl = 
(ifi) fluid/ '(7/3) sound- Quantities with over-bars are volume aver- 
aged at each iteration, and all quantities are time-averaged during 
stationary evolution. 



behavior is universal, i.e. statistically independent of 
energy injection and dissipation effects at the largest 
and smallest scales respectively. This range is referred 
to as the inertial sub-interval, and is often identified by 
self-similar behavior in the power spectrum P(k) of the 
velocity field . The 1941 theory of Kolmogorov (K41) 
(jKolmogorovl [1941) for incompressible, non-relativistic 
turbulence predicts that P(k) cx fc~ 5 / 3 in the inertial in- 
terval. In relativistic flows the amplitude of velocity fluc- 
tuations cannot exceed c, so the power spectrum of veloc- 
ity cannot obey a power law to arbitrarily large scales. 
However, the physics of a relativistic cascade may still 
be universal, and it could be revealed by an appropri- 
ate choice of generalized diagnostics, such as the power 
spectrum of fluid four-velocity. 
Figs. [2£i and[2p show the power spectrum P(k) of four- 



Spectral and Intermittency Properties of Relativistic Turbulence 



3 





Figure 2. (a) Power spectrum of the four- velocity field at resolutions 256 3 , 512 3 , 1024 3 , and 2048 3 (increasing in weight from gray to 
black), (b) The same data as in (a) but compensated by fc 5 / 3 and offset, stretched to exaggerate deviations from a 5/3 law. An arbitrary 
vertical offset is given to each curve in order to clarify the spectral shape between k/2n(= 1/L) = 10 and 100. (c) Power spectrum of the 
Hclmholtz decomposed four-velocity field. The compressive component P c (k) (black) follows a power law with index 1.80. (d) The same 
data as in (c) but compensated by fc 5 / 3 . 

is because a relativistic cascade contains a new dimen- 
sionless number, the Lorentz factor, at each scale. As 
a simple illustration, a cursory application of the K41 
dimensional argument to relativistic eddies reveals that 
the Lorentz factor at a scale I satisfies 



velocity at lattice sizes of 256 3 , 512 3 , 1024 3 , and 2048 3 , 
where 

p(k)dk = J2 fi k • K (i) 

kedfe 

is computed from the discrete Fourier transform 

< = -W £ «"(*<.m,»)e- <k - x '-" (2) 



N 3 

l,m,n 

of the four- velocity field. Throughout the text, boldface 
is used to denote the spatial components of a four-vector. 
The normalization by TV 3 , where N is the number of lat- 
tice points per direction, guarantees that P(k) satisfies 
J P(k)dk = (u-u). We find evidence for an inertial inter- 
val of relativistic velocity fluctuations between 1/10 and 
1/100 of the largest scale. As the resolution increases 
from 256 3 to 1024 3 , a short interval obeying the 5/3 law 
emerges. But the subsequent resolution of 2048 3 reveals 
a more featureful spectrum of four- velocity, consisting of 
a broken power-law which is steeper (1.92) at the largest 
scales and shallower (1.56) at moderate scales. This may 
be due to the bottleneck effect which is characterized by 
an accumula tion of power at the small-scale end of the in- 
ertial range (Falkovich 1994). There has been substantial 
effort to distinguish tru e inertial scaling from the co ntam- 
ination of bottlenecks (jBeresnvak fc Lazarian|[2009l ) . We 
report that the power spectrum of four- velocity is broadly 
5/3, owing to the fact that each higher resolution adds 
additional scales which, on the average, obey 5/3 scal- 
ing. In relativistic turbulence, asymptotically converged 
scaling behavior may not yield a single power law. This 



( 7£ ~1) 3 ( 7£ + 1) 
7? 



= Ce 2 £ 2 /c 6 



(3) 



where e is the energy injection rate per unit mass and C 
is analogous to the Kolmogorov constant. This relation 
does not admit a single power law solution. Instead, 
solutions are 7^ oc £ for large 7, and j£ oc £ 2 ^ 3 for small 
7- 

Unlike the total fluid four-velocity, the compressive 
four-velocity component it£f is strictly self-similar over 
the inertial interval. This is evidence that a local cas- 
cade of acoustic modes is operating across those scales. 
The Hclmholtz decomposition of spatial components of 
four-velocity is done in the Fourier domain, with the 
compressive part u k = Uk • k/fc and the solenoidal part 
u k = Uk — u k . As shown in Figs. [2k and [2ji, the 
power-law range of dilatational power spectrum P c (k) 
extends through the energy injection scale. This is be- 
cause only the vortical modes are being directly excited. 
Its slope is 1.80, which is in contrast with the super- 
sonic limit of highly compressible turbulence, where all 
the compressive power is in shocks, and the dilatational 
power spectrum P c (k) oc k~ 2 . The power spectrum of 
solenoidal P s (k) four- velocity is not at all self-similar, 
and accounts for the majority of power (82%) through- 
out the cascade. We remark that thermally relativis- 
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Figure 3. The generalized structure functions of four-velocity 
(longitudinal projection, squares) and Sz^t) (circles) as a 
function of the separation I for resolution 2048 3 . The vertical gray 
bar marks the relativistic scale £ 7 = 0.06L, below which velocity 
fluctuations become sub-relativistic. 

tic turbulence admits a novel mechanism for energy ex- 
change between the vortical and compressive motions of 
the fluid. The mechanical work done in compressing a 
fluid volume goes into internal energy, which for a rel- 
ativistic gas enhances its inertia through the coupling 
of total enthalpy to fluid four-velocity expressed by the 
term phu^u u in the energy-momentum tensor. This ad- 
ditional mechanism of energy transfer adds qualitatively 
new dynamics to the turbulent cascade. In particular, it 
may break the statistical decoupling between the com- 
pressive and vortical cascades that was recently d iscov- 
ered for non- relativistic compressible turbulence (|Aluiej 
1201 lh . The signature of this effect is proposed to be an 
enhanced coherency among the phases of similar-scale 
shearing and dilating velocity modes, that can be deter- 
mined from the bi-spectrum of fluid four- velocity. 

Modern descriptions of turbulence have recognized 
that two-point correlations of velocity are not normally 
distributed around the mean f l uctua t ing fluid velocit y 
(IKolmogorovl fl962l : IKraichnanl Il990t ICao et al.l fl996h . 
This effect is referred to as intermittency, and represents 
an important departure from the K41 theory. Intermit- 
tency may be quantified by examining higher moments 
of the probability distribution function (PDF) of the ve- 
locity increments at various scales £. In particular, the 
1994 model of She-Leveque (SL94) (IShe fc Leveaud 19941) 
has had great success in predicting the scaling exponents 
C p for the velocity structure functions S p (£) of order 
p. We have found that relativistic turbulence also dis- 
plays intermittency, and that the SL94 model provides an 
excellent description, provided that velocity increments 
are properly defined to accommodate relativistic fluctua- 
tions. We generalize the total velocity structure function 
as 

(4) 
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Figure 4. Dependence of £ p on the exponent p for the power-law 

scaling of longitudinal Sp (top) and total S p (£) (bottom) structure 
functions. £ p is obtained by a least-squares fit of the structure 
functions between I = 1/10 and 1/100 the domain. All the scaling 
exponents are normalized by £3. Error bars indicate the standard 
deviation over 6 equally-spaced stationary snapshots at 1024 3 . The 
intermittency models of K41 and SL94 are shown in dashed and 
solid lines respectively. 

where Su^ = u^ — u^ and £ = (Sx^Sx^) 1 ^ 2 . Four-velocity 
pairs are chosen to be simultaneous in the global center- 
of-momentum frame so that Sx° — 0. Note that S p (£) = 
(\2( lrel - l)|f/ 2 ) where "f re i is the Lorentz factor of one 
fluid element as measured in the rest frame of the other. 
S p (£) becomes (|v 2 — vi| p ) in the non-relativistic limit. 
We also generalize the longitudinal structure function 



(5) 

where £ = £/£, 



£~ to which 



S p (£) = (\Su»5u^ 2 ) 



which can also be written as (\6u ■ ■ 
because of pair simultaneity. 

We utilize to determine the scale 

relativistic velocity fluctuations persist. The kinematic 
effects of relativity are important when the relative four- 
velocity 7/3 between fluid elements exceeds 1, i.e. {^rei) > 

V2 (or S 2 {£) = 2V2 - 2 « 0.83). As shown in Fig. G2 
this occurs at £ 1 = 0.06, or 1/16 of the outermost scale. 

6*2 (£) shows two regimes of power-law scaling, having ex- 
ponents of C2 = 0.98 below £ 1 and (2 = 1-09 above £-y. 

The longitudinal structure function sj {£) obeys a power 

law with index d = 0.89 between 1/10 and 1/100 of 
the domain. This is steep relative to the K41 theory for 
incompressible non-relativistic turbulence that predicts 
(2 = 2/3, but is in close agreement with what was found 
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in simul ations of highly comp ressible non-relativistic tur- 
bulence (jKritsuk et alj|2007l ). 

We find that the SL94 intermittency model extends 
to relativistic turbulence when the structure function is 
given by Eqn. |4j Fig. |4] shows that the normalized 

exponents Cp/d through p — 6.5 are well described by 
the SL94 relation ( p = p/9 + 2 - 2(2/3) p / 3 . Error bars 
are obtained from the standard deviation over 6 snap- 
shots at 1024 3 . The total structure function given in 
Eqn. which measures the statistics of pairwise relative 
Lorentz factor, exhibits a far greater degree of intermit- 
tency. We propose that a new intermittency model for 
the statistics of relative four-velocities is required, not- 
ing that these statistics are inherently non-Gaussian from 
consideration of relativistic kinematics alone. That is to 
say, a Gaussian random velocity field does not yield a 
Gaussian PDF of j re i, owing to the non-additive nature 
of Lorentz boosts. A subsequent investigation will exam- 
ine the shape of the S p (£) PDF at fixed scales in order 
to more specifically characterize the source of the inter- 
mittency seen in Fig. |4] 



3. DISCUSSION 

In this Letter we have described a numerical investiga- 
tion of isotropic turbulence in a kinematically and ther- 
mally relativistic fluid. It is found that relativistic fluc- 
tuations persist to 1/16 of the outermost scale. Spectral 
analysis reveals an inertial sub-range of relativistic veloc- 
ity fluctuations with a broadly 5/3 index. The compres- 
sive component of the four-velocity field obeys power- 
law scaling with index 1.80 over nearly two decades. 
The SL94 phenomenology successfully describes inter- 
mittency in relativistic turbulence, provided the correct 
definition for the velocity increments is used. 

In closing, we remark that for sufficiently large 
Reynolds number, £ 7 is larger than the dissipation scale 
77 so that kinetic energy must be transmitted through a 
sub-relativistic cascade before being thermalized. In this 
case the ratio l^jL may be independent of the dissipation 
mechanism, and should thus be dictated solely by the 
Lorentz factor of the energy bearing eddies. The exact 
functional form of this dependence could be constrained 
by a parameter study over the outer scale Lorentz factor. 
Doing so would yield insight into the nature of astrophys- 



ical turbulent cascades operating far into the relativistic 
regime. 
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